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Abstract 

An analysis of the influence of the Landau-Migdal-Pomeranchuck (LPM) ef- 
fect on the development of air showers initiated by astroparticles, is presented. 
By means of computer simulations using algorithms that emulate Migdal's the- 
ory, including also the so-called dielectric suppression, we study the behaviour 
of the relevant observables in the case of ultra high energy primaries. We find 
that the LPM effect can significantly modify the development of high energy 
electromagnetic showers. 
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1 Introduction 



There are some effects that drastically reduce the cross sections of electron bremsstrahlung 
and pair production |IJ of atmospheric air showers initiated by high energy astroparti- 
cles. These suppression mechanisms can affect air showers by lengthening the showers, 
and consequently moving the position of the shower maximum deeper into the atmo- 
sphere. Two suppression mechanisms have been studied in the present work, namely, 
the Landau-Pomeranchuk-Migdal (LPM) effect [§, |3| due to the multiple scattering 
and the dielectric suppression [|lj due to the interaction of the bremsstrahlung photons 
with the atomic electrons in the medium through forward Compton scattering. 
By means of numerical simulations we have studied the influence of these effects on the 
development of air showers. We have used the AIRES program to perform such sim- 
ulations. AIRES |4[] represents a set of programs to simulate atmospheric air showers 
and to manage all the associated output data. The physical algorithms of the AIRES 
system are based on the realistic procedures of the well-known MOCCA program ||. 
AIRES provides some additional features, for example: The Earth's curvature is taken 
into account allowing safe operation of all zenith angles; the simulation programs can 
be linked to several alternative hadronic collission models; etc. To complete our study 
we have incorporated new LPM and dielectric suppression algorithms to the AIRES 
code. 



2 Migdal theory 

The LPM effect was first predicted by Landau and Pomeranchuk some 40 years ago. 
Migdal provided the corresponding quantum mechanical theory giving analytical 
expressions for the bremsstrahlung and pair production cross sections. Recently an 
experiment performed at SLAC || using targets of different compositions, measured 
the LPM effect founding that there is acceptable agreement between the experimental 
data and the Migdal's theory that is presently considered the standard treatment. 
Let us consider first the case of bremsstrahlung where an electron or positron of energy 
E emites a photon of energy k in the vecinity of a nucleus of charge Z. The Migdal 
cross section for this process is (c = 1, % — 1) 0: 

%i = ^t^GM + 2[1 + (1 - »)«]«.)>* to g*) (1) 

where 
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Figure 1. 

Bremsstrahlung 
normalized probability 
versus photon energy, as 
given by Migdal's theory, 
with and without 
dielectric suppression. 
The electron energy is 
10 18 eV and the medium 
is air (atmospheric depth 
1000 g/cm 2 ). 



= 12s 2 Qf°° e" si coth(t/2) sin(st)dt) - 6vrs 2 , (5) 

!2 if s < si 

1 + ln(s)/ln(sl) ifsi<s<l (6) 
1 if s > 1 

(si = Z 2 / 3 /184 2 ). r e is the classical electron radius (r e = e 2 /m). E LPM is the charac- 
teristic energy of the LPM effect and is given by 

_ m 4 X 

^lpm - — £2-. (7) 

Xo is the radiation length and 

£ s = m^Jin/a = 21.2 MeV. (8) 

To measure the strengh of the effect it is convenient to introduce the suppression factor 
through 

_ da LPM /dk 
da BH /dk 

where dasuldk stands for the "classical" bremsstrahlung cross section given by the 
theory of Bethe and Heitler J/J . 

The strenght of the effect largely depends on the variable s of equation (3). For 
s <C 1, the suppression is important (S <C 1), while for s 3> 1, there is no suppression 
(S = 1). In fact, when s — ► oo, the Migdal cross section reproduces, up to 3 %, the 



3 




0.3 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

u u 



Figure 2. Migdal formulation: Pair production, a) P(k,u), X = 50g/cm 2 . b) 
P(k,u), X = WOOg/cm 2 . 

main terms of the Bethe-Heitler equation Jl], [7j. 

The characteristic energy E LPM gives the energy scale where the effect is significant. 
Notice that Elpm diminishes when the density of the medium is enlarged. Therefore, 
for dilute media the LPM effect will be appreciable only for very high energies. For air 
in normal conditions, for example, taking Z = 7.3, p = 1.2 Kg/m 3 we have Elpm = 223 
PeV. 

On the other hand, when A; — ^ 0, it is necessary to take into account the change of 
the photon momemtum due the fact that the dielectric constant of the medium is 
different from one. By considering k ^> uj p one obtains: 6 = 1 — uJl/k 2 , where u p 
is the well-known plasma frequency (for air hu p = 0.73 eV). The Migdal approach 
takes into account this effect usually called dielectric suppression. The influence of 
the dielectric suppression on the bresstrahlung cross section is well noted in Figure 1 
where the normalized probability for bresstrahlung is plotted against y in the case of 
E = 10 18 eV. Notice how the emission probability is suppressed for y < 10 . Since 
the energy E = 10 18 eV this corresponds to photon energies k < 0.01 TeV. 

We go now into the pair production processes. In this case the cross section coming 
from Migdal's theory reads 

daLn,[l d E e+e ' ] = ^i G « + 2 i« + - "Oirtw* ln O < 10 > 

where 

it = 

k 



«=! (ii) 



and 



ELPAdk 
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Figure 2 shows the normalized probability of pair production at atmospheric depth of 
50 and 1000 g/cm 2 and for different energies of the primary photon. Notice how the 
production probabilities are progresively suppressed when the primary energy rises. 
From the figure above, it is also evident that symmetric processes are more affectted 
by the suppression mentioned. 

3 Practical Implementation 

The AIRES simulation system @] has been used as a realistic program to perform the 
simulations needed to make our analysis. The LPM effect and the dielectric suppression 
have been incorporated into the AIRES program. 

When the characteristics of the atmosphere are taken into account it comes out that 
the LPM effect must be considered for all the events where the energy of the primary is 
larger than 100 TeV. This ensures that the effective cross sections are calculated with 
a relative error that is never larger than a few percent. 

4 Simulations 

We have analized the influence of the LPM effect and the dielectric suppression on the 
air shower development, performing simulations for different primary energies (from 
10 14 to 10 21 eV) and different primary particles initiating the showers (gamma, proton, 
electron) . 

For showers initiated by gammas with primary energy larger than 10 20 the impact 
of the LPM effect is evident. It affects the position of X max (the maximum of the 
shower), the number of particles at X max and the fluctuations of these magnitudes. 
This is shown in Figure 3, where the total number of charged particles is plotted 
against the vertical depth (longitudinal development of all charged particles) for two 
different primary energies. The LPM effect affects the gamma showers lengthening and 
consequently moving the average position of X max deeper into the atmosphere. This 
effect is evident in Figure 4. a where one sees that X max is shifted in approximately 100 
g/cm 2 for 10 20 eV and 500 g/cm 2 for 10 21 eV. The fluctuactions of A max also increase 
when the LPM effect is taken into account as can be seen in Figure 4.b. The average 
number of charged particles at X max is reduced and its fluctuations are larger when the 
LPM effect is introduced as can be seen in Figure 5 for gamma showers of more than 
10 20 eV. We have not found appreciable differences for gamma initiated showers with 
primary energies less than 10 18 eV. This allow us to conclude that even if the LPM 
effect must be taken into account for all particles with energies larger than 100 TeV, 
the fraction of such events is statistically significant only for electromagnetic showers 
initiated by primaries with energies larger than 10 18 eV. 
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In Figure 6 the total number of ^ is plotted against the vertical depth (longitudinal 
development of /jl ) againts two different energies. This observable is very important 
in the determination of the primary composition. In agreement with the longitudinal 
development of all charged particles, it presents the same modifications when the LPM 
effect is taken into account. 
For air showers initiated by protons, no measurable differences appear neither in the 
average X max , nor in the fluctuations of this quantity. The protons interact hadroni- 
cally, and then, the electromagnetic shower, where the LPM effect takes place, starts 
later. The proton primary energy is shared among the secondary showers after the first 
interaction and the electromagnetic cascade begins with energies that are about 2-4 
orders of magnitude less than the inicial proton energy. Therefore, one should compare 
the 10 20 proton showers with the gamma showers of initial energies of 10 17 -10 18 eV, 
where we have not found appreciable differences between showers with and without 
LPM effect. The longitudinal development of all charged particles for proton initiated 
showers is plotted in Figure 7. The differences between the simulations with and with- 
out dielectric suppression are less important and are shown in Figure 3 and 6 for the 
showers initiated by gammas. 

Finally notice that the characteristics of electron initiated showers are very similar to 
those corresponding to gamma showers and for this reason we have not included here 
any related plots. 
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Figure 3. Longitudinal development: All charged particles. Parameters: Primary 
particle: 7. Primary energy: 10 19 eV, 3 x 10 20 eV. Zenith 60°. Thinning energy: 10 
rel. Injection altitude: 100 km. 
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Figure 4. Shower maximum X„ 
proton showers. 



(a), and X max fluctuations (b) for gamma and 
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Figure 5. Number of charged particles at X max , N max (a), and N max fluctuations (b) 
for gamma showers. 
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Figure 6. Longitudinal development of muons. Parameters: Primary particle: 7. 
Primary energy: 10 19 eV, 3 x 10 20 eV. Zenith 60°. Thinning energy: 10 -5 rel. 
Injection altitude: 100 km. 
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Figure 7. Longitudinal development: All charged particles. Parameters: Primary 
particle: p. Primary energy: 10 19 eV, 3 x 10 20 eV. Zenith 60°. Thinning energy: 10 -5 
rel. Injection altitude: 100 km. 
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5 Conclusions 



The LPM effect introduces modifications on the development of gamma and electron air 
showers if the primary energies are larger than 10 19 eV. These effects can be observed 
in the longitudinal development of the showers. The X max position for such initial 
conditions moves deeper into the atmosphere and its fluctuations are increased when 
the LPM effect is taken into account. The longitudinal development of ^ also changes 
in concordance with all charged particles case. 

We have not found any significant effect if the showers are initiated by proton with 
primary energies up to 10 21 eV, because in this case, the electromagnetic shower, where 
the LPM effect takes place, begins later, when the initial energy is shared among the 
secondary particles, reducing the initial proton energy in 2-4 orders of magnitude. 
Clearly the same reasoning is valid for nuclei primary cosmic rays. 
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